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, We consider concentration dependence of rheological properties of associative telechelic polymer 

£Nj ■ solutions. Experimental results for model telechelic polymer solutions show rather strong concentra- 

tion dependence of rheological properties. For solutions with relatively high concentrations, linear 
viscoelasticity deviates from the single Maxwell behavior. The concentration dependence of charac- 
M-H teristic relaxation time and moduli is different in high and low concentration cases. These results 

suggest that there are two different concentration regimes. We expect that densely connected (well 
percolated) networks are formed in high-concentration solutions, whereas sparsely connected (weakly 
percolated) networks are formed in low-concentration solutions. We propose single chain type tran- 
sient network models to explain experimental results. Our models incorporate the spatial correlation 
effect of micellar cores and average number of elastically active chains per micellar core (the network 
functionality). Our models can reproduce non-single Maxwellian relaxation and nonlinear rheolog- 
ical behavior such as the shear thickening and thinning. They are qualitatively consistent with 
^ 1 experimental results. In our models, the linear rheological behavior is mainly attributable to the 

q , difference of network structures (functionalities). The nonlinear rheological behavior is attributable 

to the nonlinear flow rate dependence of the spatial correlation of micellar core positions. 
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Telechelic associative polymers which have hydrophilic main chains and hydrophobic chain ends form self-assembled 
micellar structures in aqueous solutions [J-Q- If the concentration of telechelic polymers is sufficiently high, the micellar 
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" cores are bridged by polymers to form a percolated, network like structure. Such a network structure shows various 
kinds of interesting rheological behavior and rheological properties of telechelic associative polymer solutions have 
been extensively studied. For example, various experimental results have been reported for hydrophobically modified 
cthoxylated urethanc (HEUR) aqueous solutions [2H5|. 
■ One characteristic rheological property of telechelic polymer solutions is the linear viscoelasticity being well de- 
$_i 1 scribed by a single Maxwellian form in a rather wide range of experimental parameters. Theoretically, such a rheo- 
logical property can be explained in the framework of the transient network model 0,0]. The characteristic relaxation 
time is considered to be the dissociation time of a chain end, and the characteristic modulus is considered to be the 
elastic modulus of the network formed by telechelic polymers. However, the simple description by the transient net- 
work model is not always valid. For relatively high concentration solutions, deviation from the single Maxwellian form 
has been report ed[8l4lO|. In most of the transient network type models, such non-single Maxwellian relaxation cannot 
be explained straightforwardly. (This is due to the nature of the mean field type approximations in the transient 
network type models.) This implies that we need to incorporate some additional factors, which arc affected by the 
polymer concentration, to the models. 

Another characteristic rheological property is that telechelic polymer solutions often exhibit shear thickening behav- 
ior. The transient network model has been improved to reproduce such experimentally observed nonlinear rheological 
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behavior. In sonic improved versions of the transient network model[3j, Illl4l7j |. the nonlinear elasticity (or finite 
extensively) and the stretch dependent dissociation rate are introduced. It has been reported that the combination 
of the nonlinear elasticity and the stretch dependent dissociation rate can reproduce several nonlinear rhcological 
properties. For example, models by Tripathi et al.\^ and by Indci, Koga and Tanaka [l4l - fl7j can reproduce several 
nonlinear features such as the shear thickening and thinning, by tuning some fitting parameters (which represent the 
strengths of the nonlinear elasticity and the stretch dependent dissociation). 

However, a recent study showed that conventionally proposed transient network models cannot explain some rhe- 
ological data of a HEUR aqueous solution[lj|, such as the first normal stress coefficient data that exhibit no 
nonlinearity in the shear-thickening regime of the shear viscosity (77). To explain this experimental fact, a new tran- 
sient network model with an anisotropic (and shear rate dependent) bridge formation process was proposed. The 
anisotropic bridge formation model can qualitatively explain this behavior of as well as the shear thickening (and 
thinning) behavior of rj under fast shear, the latter being explained also by conventionally proposed models. Although 
the new model provides a new viewpoint to understand dynamical behavior of telechelic polymer solutions, it was 
constructed in a phenomenological way and the molecular level mechanism of anisotropic bridge formation has not 
been specified explicitly. 

From the viewpoint of nonequilibrium dynamics of flowing systems, anisotropics under flow are physically rea- 
sonable. For example, colloid dispersions show structural anisotropies and nonlinear rhcological behavior under 
flow [l9l I20I ]. The interaction potential and spatial correlation between colloid particles are important to under- 
stand/analyze such systems. In telechelic polymer solutions, there are also interactions and spatial correlations 
between micellar cores || [Io|, HH, [H[ (although their forms are quite different from colloid dispersions). Clearly, such 
interactions or spatial correlations depend on the polymer concentration, and thus we expect the rhcological proper- 
ties of telechelic polymer solutions to be rather strongly dependent on the concentration. Indeed, we can find rather 
strong concentration dependence of several rhcological properties in literature 0, 

This implies that we should design a model for telechelic polymer solutions which directly takes into account the 
spatial correlation (or the interaction potential) of the micellar cores, to naturally reproduce concentration dependence 
of rheological properties and anisotropic correlations and nonlinear rheological behavior under shear. Nevertheless, 
the importance of the spatial correlations for the rheology of telechelic polymer solutions has been rarely pointed out. 
Especially, the spatial correlations have been overlooked or ignored for single chain type models, as far as the authors 
know. 

In this work, we conducted linear viscoclastic measurements for HEUR aqueous solutions with different concentra- 
tions. We show that the linear rheological behavior is strongly affected by the concentration. The experimental data 
suggest that there are two different concentration regimes. Then we propose a simple but non-trivial single chain type 
transient network models for telechelic associative polymers, which explicitly take into account the spatial correlation 
information. We show that our models give both single and non-single Maxwellian relaxation moduli depending on 
the polymer concentration. We also show that nonlinear rheological behavior such as the shear thickening and thin- 
ning can be reproduced by our models considering neither nonlinear elasticity nor stretch dependence of the bridge 
reconstruction rates. We discuss concentration dependence of rheological data from the view point of the transient 
network type models. 

II. EXPERIMENTAL 
A. Materials and Measurements 

HEUR with hcxadccyl end groups and polyethylene oxide (PEO) main chain was used. The HEUR was synthesized 
previously [l8[ with a conventional method, and the molecular weight and polydispersity index were M w = 4.6 x I0 4 
and M w /M„ = 1.35 (the PEO precursor molecular weight and polydispersity index were M w = 1.9 x 10 4 and 
M w /M n = 1.1). The details of the synthesis and characteristics are found in Ref. (l8j | . 

The HEUR aqueous solutions were prepared by dissolving HEUR into distilled water through stirring for 24 hours. 
For comparison, the PEO aqueous solutions with a similar molecular weight and the same concentrations as the 
HEUR solutions were also prepared. The PEO sample was purchased from Aldrich (M w = 3.9 x 10 4 , M w /M n = 1.08) 
and dissolved into distilled water. 

Storage and loss moduli, G'(ui) and G"(oj), were measured for HEUR aqueous solutions at various temperatures 
and concentrations. The measurements were conducted with a stress-controlled rheometer (MCR-301, Anton Paar). 
A cone-plate type fixture with the diameter d = 75 mm and the cone angle 9=1° was utilized. The measurement 
temperatures were T = 5, 10, 15, 20, and 25°C, and the concentration c was varied from lwt% to 10wt%. The angular 
frequency range was set as 0.03 rad/s < tu < 500 rad/s. Storage and loss moduli of PEO aqueous solutions with the 
same concentrations and the similar molecular weight as HEUR solutions were also measured in the same condition. 
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(Due to the low moduli of solutions and the limitation of the equipment, only the loss moduli data in a limited range 
1 rad/s < u < 10 rad/s were obtained with acceptable accuracy.) 



B. Results 



First we show the master curves for the storage and loss moduli, G'(uj) and G"(ui) of lwt% HEUR aqueous 
solution at the reference temperature T r = 25°C in Figure [TJ The master curves (master plots of b^G' and b^G" 
against ar,ooj) are constructed by shifting G'(w) and G" {ui) data at different temperatures horizontally and vertically, 
to superimpose the low frequency data. The horizontal and vertical shirt factors are denoted by ot.o and br,o, 
respectively. As reported in the previous work[l8j, the time-temperature superposition works well for this HEUR 
solution and the G'(u>) and G"(ui) data are fitted well by a single Maxwell model, except the high frequency region. 

G '<"> = G °TT^r g "<"> = g "T^F (1) 

Here r and Go are the characteristic relaxation time and modulus, respectively. The deviation in the high frequency 
region means that there is another relaxation mode of which temperature dependence is different from the relaxation 
mode in the low frequency region. In the followings, we call these modes in the high and low frequency regions as 
fast and slow modes, respectively. 

Ng et aZ.[23[ proposed that the fast mode is due to the relaxation of micellar cores, whereas Bedrov et a/.[24| claimed 
that the fast mode corresponds to the local relaxation of hydrophilic main chains. If the fast mode corresponds to the 
local relaxation mode of PEO chains, it should show the same temperature dependence as a PEO aqueous solution 
without associative end groups. We determined the horizontal shift factor a^.i of the lwt% PEO solution (utilizing 
the vertical shift factor &t,i = T/T r ). Then we shifted the G' and G" data of the lwt% HEUR solution with the 
shift factors CbT,i and &t,i- The result is shown in Figure [5J and ar,i is compared with a^.o 01 the HEUR solution 
in Figure [3] As shown in Figure [U ax,i allows the G'(w) and G"(w) data to collapse into master curves at high 
frequencies where the fast mode is observed. a<r,i is less dependent on T compared with axfi. This result means that 
the fast mode corresponds to the local relaxation of PEO chains, as claimed by Bedrov et aZ.(24|. However, it should 
be also noticed that the relaxation time and the characteristic modulus of the fast mode of the HEUR solution are 
much larger than those of the PEO solution with the same concentration, mostly because the HEUR chains form 
an associated network-like structure and the relaxation modes of such a network are different from ones of the free 
chains. We will discuss this point later. Although we have shown the master curves only for the lwt% solution case, 
we can similarly construct master curves for slow and fast modes for other concentration cases. 

We show the storage and loss moduli data at T = 25°C for different concentrations (c = 1,2,5, and 10wt%) in 
Figure |U We can observe slow and fast modes for all concentrations. However, the w dependence of G'(w) and 
G"(w) changes with the concentration rather largely. This change is more prominent for the fast mode than for the 
slow mode. For c = 1 and 2wt%, the G"(w) at high uj (the fast mode) is almost proportional to u>, which indicates 
that the fast and slow modes are rather separated. For c = 5 and 10wt%, on the other hand, the G"{uj) at high u> 
becomes less dependent on w, and the fast and slow modes are not well separated compared with the cases of the 
1 and 2wt% solutions. The slow mode distribution also changes with the concentration, although it is not so clear 
in Figure |31 To see this change clearly, we show the G'(w) and G"(w) rescaled by the characteristic time r and 
modulus Go for the slow mode in Figure [SJ It is evident from Figure [5] that the slow mode deviates from the single 
Maxwellian form as the concentration increases. Similar non-single Maxwellian relaxation moduli have already been 
reported in some experimental or simulation studies pl-liol [25Tj29j . For example, Berret and coworkers (8l-[Toj reported 
that the shear relaxation moduli of solutions of F-HEUR (HEUR with pcrfluoroalkyl chain ends) can be fitted by the 
stretched exponential form, G(t) = Go exp[— (t/r) 01 ] (with Go, r and a being the characteristic modulus, characteristic 
relaxation time, and the power- law exponent). Mistry et al. reported non-single Maxwell relaxation modulus data for 
relatively high concentration poly(butylene oxide)- (etylene oxide)- (butylene oxide) triblock copolymer solutions |30j. 
Some shear relaxation modulus data for relatively high concentration HEUR solutions 0, H| also deviate from the 
single Maxwellian form. 

Finally we show the characteristic relaxation time r and modulus Go data of the slow mode in Figure [6] Both r and 
Go depend on the HEUR concentration. The concentration dependence at low and high concentration regions can be 
individually fitted to power laws. The characteristic time r depends on the concentration casra c 0-63 anc j T ^ c o.32 
for low and high concentration regions, respectively (Figure (Ha)). Similarly, the characteristic modulus depends on c 
as Go oc c 2 ' 3 and Go oc c 18 (Figure |6|b)). Both r and Go data show stronger c-dependence at the low concentration 
region. From Figure HI the crossover concentration c c is roughly estimated as c c « 4wt%. At c > c c , the G'(w) and 
G"(co) data deviate considerably from the single Maxwellian forms, as noted in Figure [SJ Here it should be noted 
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that Annable et al. [2J already measured the concentration dependence of r and Go of HEUR solutions systematically, 
and reported similar data. (They proposed a statistical model to explain their experimental data. Their model takes 
account of superbridge and superloop structures, and can reasonably reproduce the concentration dependence of r 
and Go.) Further discussions of the experimental results in relation to conventional transient network models will be 
made in Sec. [IVJ 



III. THEORETICAL MODEL 



As we have shown in Section UH the linear viscoelasticity of a telechelic polymer solution depends on the concentra- 
tion and there are two different concentration regimes. The nonlinear rhcological data of the lwt% aqueous solution 
of the same HEUR has suggested the existence of the anisotropic bridge formation dynamics under flow fig . To 
explain such rheological behavior by theoretical models, we need to employ models which can take into account the 
polymer concentration and some spatial correlation between micellar cores. One possible candidate is the statistical 
model by Annable et al.fy, which can reproduce concentration dependence of linear rheological properties. However, 
it is not easy to study nonlinear rheological properties by their model. We consider that the transient network type 
models are suitable to study both the linear and nonlinear rhcological properties. In this section we propose transient 
network type models which take into account effects of the average number of bridge chains per micellar core and the 
correlation between micellar cores. We show that these factors qualitatively affect several rheological properties. We 
consider only the slow mode, and thus the models in the followings can not reproduce the fast mode. 



A. Dense Network Model 



Wc first consider cases where the average number of bridge chains per micellar core (the functionality /) is sufficiently 
large (typically / > 3) and dense networks are formed. A schematic image of a well-percolated network is shown in 
Figure [7|a). We assume that the polymer chains are not sufficiently long and not too much concentrated so that the 
entanglement effects are negligible. For such cases, it will be sufficient to consider only one tagged (labeled) chain 
in the system. In the followings, we call such networks as "the dense networks" and the model described in this 
subsection as "the dense network model" . 

Before we consider the dynamics, we focus on the equilibrium probability distribution function for a single, tagged 
chain in the system. We assume that a polymer chain can take two states; the bridge and the loop. (The dangling 
chains are not explicitly considered here because the characteristic rcassociation time scale of the dangling chains 
should be much shorter than those of the bridge or loop states.) A schematic image is depicted in Figure EJh). Wc 
express the state of the chain by a state variable n, which takes (loop) or 1 (bridge). If the chain is in the bridge 
state, the information of the end-to-end vector r is needed to fully specify the state. We can interpret r as the relative 
position of the partner core, as shown in Figure B^b). 

Under the mean field type approximation, the interaction between two micellar cores is expressed by an effective 
mean field interaction potential v(r). The effective interaction potential consists of several different effects such as 
the steric repulsion (due to the corona chains) and hard core like interaction between micellar cores. (We expect that 
the steric repulsion will be dominant, and practically other contributions may be ignored.) If a bridge is formed, the 
chain feels an extra potential energy. Roughly speaking, this extra energy corresponds to the elastic free energy of a 
chain, u(r), and the equilibrium probability distribution is given as 



f 



P eq (0) = -Z 

Qi 

p cq(l,r) = exp 



v(r) + u(r) 



knT 



(2) 
(3) 



Here A is the thermal de Broglie wave length, and S is the partition function. Z and Q\ are the partition function 
of a single loop chain and the partial partition function of a single bridge chain, respectively. The quantities with the 
subscript "eq" represent the equilibrium quantities. 



S = 2o 



|f / dr ex P 



v(r) + u(r) 



(4) 



The effective interaction potential can be related to the structure of the multi chain system. It should be determined 
so that the resulting spatial structure becomes consistent with the structure of a target system, such as the structure 
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detected by scattering experiments. Under the mean field approximation, the radial distribution function (RDF) of 
the system can be easily calculated. 



Po g cq (r) = n P cq (l,r) + p e- v( - r ^ kBT P cq (0) 



PoZo,i , s.-u(r)/k B Ti -v(r)/k B T ^ 



Here po is the average spatial density of the micellar core, no is the aggregation number of a miccllar core, and 
£ = noQi/ pqA 3 Zq is the effective activity (fugacity). We have assumed that the system size is sufficiently large. Eq 
([5]) gives the relation between the mean field potential and the radial distribution function. 

v(r) = -k B T\ng cq {r) + fc B Tln[l + ^- u{r ^ kBT ] (6) 

Here we have dropped an additive constant which does not affect thermodynamic properties. The first term in the 
right hand side of © is the effective potential form that appears under some approximations (for example, in the 
liquid state theory [3 If). The second term in the right hand side of ((6]) represents the repulsive interaction between 
cores which cancels the attractive interaction by bridges. We cannot reproduce the correct spatial correlation without 
this repulsive interaction. 

For convenience, we introduce another equilibrium probability distribution function. We express the (unnormalized) 
conditional probability distribution of a micellar core under the condition n = as 

$ oq ( r ) = Pl e -v(r)/k B T (7) 

rio 

It is straightforward to show that the following relation holds. 

P cq (l, r) = Ce-«W/ fc « r $ eq (r)P eq (0) (8) 

Eq |7]) is useful when we consider the bridge construction process. This is because $> eq (r) represents the probability 
to find a partner core at a certain position r in space, which is required when we consider the bridge construction 
process. 

To study dynamical properties such as rheological properties, we need the dynamic equations for probability distri- 
butions. We consider the time evolution equations for three time-dependent distribution functions (P(l, r, t), P(0, t), 
and $(r, t)). Roughly speaking, there are two different contributions for dynamic equations. One is the reconstruction 
of bridge chains, and another is the motion of micellar cores. 

First we consider the bridge reconstruction dynamics. So far, several different destruction rate models have been 
proposed for transient network type models pit Pzl. [l3l [l4| . Here we describe the destruction rate in the following form. 

W(0|l,r) = -L (9) 
r(r) 

r(r) represents the characteristic destruction time for the bridge chain with the end-to-end vector r. (In the next 
section, we simply set r(r) as a constant.) The bridge formation rate should be proportional to the probability that 
we find a partner core. Then the bridge formation rate can be modeled as W(l, r\0)§(r, t). From the detailed balance 
condition, the explicit form of W(l, r|0) is automatically determined. 

In most of conventional transient network models, the dangling-to-bridge transition rate is utilized instead of the 
bridge construction (loop-to-bridge transition) rate. This difference is not serious, because the dangling-to-bridgc 
transition rate can be cast into the bridge construction rate, unless the system is subjected to a very fast flow. 

Here it is worth noting that in many transient network models, both the dangling-to-bridge and bridge-to-dangling 
transition rate models have been proposed and utilized to reproduce rheological data well. In other words, both two 
transition rates are freely tunable in these models. However, as long as the equilibrium probability distribution is given, 
the bridge construction rate is automatically determined and no longer freely tunable. (If the bridge construction rate 
is given, the bridge destruction rate is automatically determined from the detailed balance condition.) Our bridge 
reconstruction rates shown above do not affect the equilibrium probability distribution. 

Next we consider the contribution of the motion of micellar cores. Micellar cores move by the thermal noise (the 
Brownian motion) and by externally imposed flows. These effects can be modeled by the (generalized) Fokkcr-Planck 
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operator. The Fokkcr-Planck operator can depend on whether the chain is in the bridge state or in the loop state, 
and thus we employ an n-dependent Fokker-Planck operator, C(n,t) (n = and 1, for loop and bridge). In or near 
equilibrium, under the Markov approximation, the Fokker-Planck operator C(n, t) can be expressed simply as 



d[nu{r) + v(r)\ p + ^ T &P(r) 



Or 



dr 



l-.[ K (t)-rP{r)\ 



(11) 



where (q is the effective friction coefficient and n(t) is the velocity gradient tensor. The effective friction coefficient 
Co can be related to the long time self diffusion coefficient of a micellar core Dq as 



Dn = 



2k B T 
Co 



(12) 



The numerical factor 2 comes from a fact that r is the distance between two micellar cores and they move almost 
independently in the long time limit. 

It should be noticed that eq (fTTj) is not always valid. For example, if the memory effect [32| is not negligible, 
we need to introduce the memory kernel. In noncquilibrium states, such as under shear flow, the dynamics will 
be drastically affected by external driving forces due to flow. If the system largely deviates from equilibrium, the 
effective potential can be different from the equilibrium potential. Also, the mobility tensor (friction coefficient) can 
become qualitatively different from the equilibrium form (33l. |3~4|. Therefore, for noncquilibrium cases, we may need 
to employ a noncquilibrium effective free energy and an anisotropic mobility tensor. It is not a simple task to model 
nonequilibrium dynamics, and in the current work we mainly considers dynamics near equilibrium. Fortunately, the 
explicit form of C(n, t) does not severely affect the results of the following analyses (at least qualitatively), and we do 
not further discuss these effects here. 

Finally, by using the bridge reconstruction rates and the Fokkcr-Planck operators, we have the following dynamic 
equations for probability distributions. 



dP(l,r,t) 

dt 
dP(0,t) 

dt ~ 
d$(r,t) 

dt~ 



C(l,t)P(l,r,t) + W(l,r\0)<f>(r,t)P(0,t) ~ W(0\l,r)P(l,r,t) 
dr [W{0\1, r)P(l, r, t) - W(l, r|0)$(r, t)P(0, t)] 
= £(0,t)*(r,t) 



(13) 
(14) 
(15) 



The dynamic equation for the probability distribution of the micellar core, &(r, t) (eq (fl~5j) ). was not introduced in the 
previous transient network models (3l Illl4l7j . (The previous models implicitly assumed that <fr(r,t) does not depend 
on flow history.) As we will show in the next section, $(r,i) largely affects rheological properties. 



B. Sparse Network Model 



If the functionality is relatively small (/ < 3), the dense network model introduced in the previous subsection is 
not appropriate. In such cases, many (though not all) micellar cores in a network link only two bridge chains and 
can not be regarded as active nodes sustaining the network elasticity. Then networks are mainly formed by so-called 
superbridge structures 0, and the (apparent) fraction of elastically active chains become small. We call such networks 
as "the sparse networks" and the model described in this subsection as "the sparse network model". A schematic 
image is depicted in Figure E^a). 

In the sparse network model, most of chains in the system is not elastically independent (See Figure E{b)), and thus 
we should use a superbridge as the elementary unit (effective bond) to construct a mean field model. For simplicity, 
we assume that the number of bridge chains which form a single superbridge is constant and express it as m {m > 2) . 
We assume the effective interaction potentials are the same as ones in the dense network model. Even with these 
assumptions, it is still difficult to accurately formulate a mean field model for a single superbridge. In the current 
work, we attempt to construct a model which captures the essential properties, rather than a quantitatively accurate 
model. As we will show later, many features of the sparse network model are consistent with experiments even if we 
employ rough approximations and simplifications. 

We express the probability that a superbridge (with the end-to-end vector r) is present as P oq (l,r). As a rough 
estimate, this can be approximately expressed in terms of the equilibrium probability distribution of a bridge chain 
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(in the dense network model, eq (J5J). In the followings, we utilize the short-hand notations for convolutions. 

[f*g](r) = J dr'f(r-r')g(r') (16) 

/rn 
dr ... dr m 5(r - r m )S(r Q ) J[ f(r k - r fc _i) (17) 



We approximately express P eq (l,r) as 



P eq (l,r)aPM(l,r)cx 



A 3r; 



k=l 



-{u+v)/k B T 



(r) 



(18) 



To make the expressions similar to the dense network model, we introduce the following effective bond potential 
u(r). 



u(r) = —v(r) — /cgTln 



1 



A.3(m-1) 



,-(u+v)/k B T 



(r) 



(19) 



Substituting eq (|19[) into cq (|18p . we can express the equilibrium probability distributions in the forms similar to eqs 
© and ©. 



P oq (0) » $ 
Aq(l,r) 



Qi 



A 3 S 



s CX P 



u(r) + v(r) 



(20) 
(21) 



where iE>o is the partition function of an inactive bond (a loop, a dangling end, or a supcrloop), and Qi is the partial 
partition function of an active superbridge. S is the partition function and defined as follows. 



u(r) + v(r) 



(22) 



Unlike the dense network model, the superbridge becomes elastically inactive when one bridge chain in a superbridge 
is destructed. The destruction rate is expected to be hardly dependent on the chain stretch (unless the flow is very 
fast), and thus we assume that the characteristic destruction time is constant. (According to the Kramers theory [32|. 
the destruction rate is determined by the curvatures of the potential at the local minima and maxima, and the energy 
barrier. Both of them are not largely affected by the chain stretch.) Thus the destruction rate of the superbridge 
becomes 



- , , , m 
W(0\l,r) w — 

TO 



(23) 



with To being the characteristic (or average) destruction time of a single bridge chain. The superbridge construction 
rate is not so simple compared with the dense network model. Because the destruction event of an superbridge occurs 
when one of the bridge chains in the superbridge is destructed, the construction occurs when (m — 1) bridge chains are 
present and a new bridge chain is constructed. Such a consideration suggests the construction rate in the following 
form. 



W'{l,r\Q) 



A 3 ( m ~ 



-u/k B Tj£\[m-k] 



$ * (e 



-«/*b t $)[M 



k=l 



(r) 



W'(l,r\0) is determined from the detailed balance condition as 



^'(l.rlO) = —{= e ~ a{r)/kB7 'elvM-vir^/ksT 



ro 



with 



v(r) = -fcsTln 
Z n^ 



A 3 < 



rn 

1 _ ^2^ c -(y.+v)/k B Ty m -k] ^ e -v/k B T ^ ( e -(u+v)/k B Tj[k--L]^ 



fc=l 



(24) 

(25) 

(26) 
(27) 
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In the fallowings, we express the summation of convolutions in eq ([24)) shortly as follows. 

. [v(r)-v(r)]/k B T m 

$H (r) EE A3(m _ 1} ^[( e -«/*BT $ j[m-fc] # $ ^ ( e -«/*BT $ j[fc-l]^ p j ( 28) 

fe=l 

Also, it is convenient to introduce the following construction rate. 

W(l, r|0) = r|0)e-[ s M- l 'Ml/ feBT = —^ e - u ^/ kBT (29) 

TO 

The construction rate is then rewritten as W(l, r|0)<f>[ m l (r), which is formally similar to that in the dense network 
model. 

Under the Markov approximation, the Fokkcr-Planck operator for a superbridge will be expressed as follows. 



d[u(r) + v(r)} dP(r) 
P{r) + KbT- 



dr dr 



A ■ ■ rP(r)] (30) 



Here Co is the effective friction coefficient for the end-to-end vector of a superbridge. This effective friction coefficient is 
related to the diffusion coefficient of clastically active nodes (to which three or more active superbridges are attached) . 
We expect that £o is larger than £ , and thus the effect of the Brownian motion will be small compared with the case of 
the dense network model. Also, from eqs (|2"3"|) and ([2"§]) . the characteristic time scale of the superbridge reconstruction 
process is relatively small. Therefore, we expect that the effect of the Brownian motion will be negligibly small for 
practical cases. 

The dynamic equations for the sparse network model are finally described as 

dP(l,r,t) = £( 1)t )p( 1)r) ^ + ^( 1)r | )g[^]( r) ^p( 0)t )_^( | 1)r )p( 1)r)t ) (31) 



at 

dP(0,t) 



dr 



W{0\l,r)P(l,r,t) - W(l,r\0)<f>l m ](r,t)P(0,t) 



(32) 



The dynamic equation of 4>(r, t) is the same as the dense network model, eq (|15|) . While eqs (|3Tj) and |32|) arc similar 
to eqs (H2J) and tH}, the dependence of the construction rate on <&(r,i) is qualitatively different: $(r,i) in the dense 

network model is replaced by a convolution $[ m l(r,t). 

Before we proceed to analyses of rhcological properties, we shortly comment on the relations among the dense and 
sparse network models, and the anisotropic bridge formation model. Although we have assumed that a superbridge 
in the sparse network model consists of multiple bridge chains (m > 2), it is formally possible to set m = 1. In the 

case of m = 1, convolutions simply become /' m '(r) = f(r) and $[ m l(r) = $(r). For this case, the sparse network 
model reduces to the dense network model with r(r) = r . 

However, for m > 1, the sparse network model is qualitatively different from the dense network model. We need to 
employ the sparse network model if the functionality is small and the contribution of superbridges is not negligible. We 
expect that the dependence of the superbridge construction rate on the spatial correlation is enhanced in the sparse 
network model. This would be intuitively natural because in the sparse network model, one superbridge consists of 
multiple bridge chains connected in series and the superbridge conformation strongly reflects the anisotropy of the 
underlying micellar core distribution. 

The sparse network model reduces to the previously proposed anisotropic bridge formation model, under some 
conditions. (The details are shown in Appendix [AT) If we compare dynamic equations (f3"Tj) , (|3"2"|) and (TT5)) with the 
anisotropic bridge formation model[l8[, we find that $(r, t)P(0, t) works as the source function for a newly constructed 
bridge. The source function in the anisotropic bridge formation model was introduced rather phcnomcnologically and 
its molecular meaning was left rather arbitrary. In our present work, the source function is directly related to the 
spatial distribution of the micellar cores. Therefore, models formulated here lend support from a molecular point of 
view, to some extent, to the use of the phcnomcnologically designed model. 



C. Rheological Properties 

In this subsection we show some rheological properties calculated from our models shown in Sections IIII Al and 
IIIIBI To make the models simpler and tractable, in the followings we limit ourselves to simple systems. We assume 
that effects of the nonlinear elasticity and the stretch dependent bridge reconstruction are absent. Although these 
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effects arc considered to be essentiai in some transient network models, they are not always necessary to reproduce 
experimentally observed rheological properties. (We will discuss the effects of nonlinear elasticity and the stretch 
dependent bridge reconstruction rates later.) 

We employ the constant bridge destruction time model 



r(r) = r 



(33) 



where to is a constant independent of r. (This assumption was already introduced for the sparse network model.) We 
also employ the harmonic elastic potential of a bridge chain (the Gaussian chain model) . 



u(r) 



ik B T 



(34) 



where Rq is a constant which corresponds to the average end-to-end distance of a chain. 

To calculate rheological properties, we need the microscopic expression of the stress tensor. In this work we assume 
that the stress is mainly caused by the bridge chains. For the dense network model, we define the microscopic stress 
tensor as 



<x(n) 



n p 







du(r) 
dr 



r - k B Tl 



(n = 1) 
(« = 0) 



(35) 



Here, uqPq/2 is the number density of polymer chains and 1 is the unit tensor. This expression is consistent with the 
stress-optical rule because we have employed a harmonic potential for u{r) (eq (|34|) ). It should be noted that eq (|35|) 
is not a unique candidate. For example, we may define the stress tensor in a way that it becomes a conjugate variable 
to the deformation (In this case, we have an extra term which is proportional to [dv(r)/dr]r. However, this term is 
not important unless the micellar cores are rather concentrated.) For a given set of probability distribution functions 
P(l,r), -P(O), and $(r), the average stress tensor is calculated as (cf. eqs (fM]) and (pJ5j) ) 



tr = J dr&(l)P(l,r) 



n Q p Q k B T 



j dr 


3rr 






. R l 



P(l,r) 



(36) 



For the sparse network model, we define the stress tensor as follows, instead of eq (l35|) . 

du (r) 



n c gp 




Or 



-r - k B Tl 



(n = 1) 
(n = 0) 



(37) 



Here n e ff is the effective number of elastically active supcrbridges per micellar core, and n G ffpo represents the effective 
number density of supcrbridges. The relation between n e g and n is generally not simple, but physically n c ^ should be a 
monotonically increasing function of the polymer concentration and it should approach no/2 in the high concentration 
limit. 

Although there are qualitative differences between the dense and sparse network models, their equilibrium probabil- 
ity distributions or dynamic equations are apparently similar. Thus, in the followings, we mainly describe the detailed 
derivations and analyses for the dense network model, which is simpler than the sparse network model. The results 
for sparse network model are obtained in a similar way, and we just show the results without detailed derivations. 



1. Shear Relaxation Modulus 



Here we calculate the shear relaxation modulus G(t) of the dense network model. We assume that the system is 
in equilibrium at t < 0, and a small step shear deformation is applied at t — 0. For convenience, we rewrite the 
Fokkcr-Planck operators (in eq (JTTJ) ) as 

C(n, t)P(r) = £„(n)P(r) - A • fr5(t)r y P(r)} (38) 

where £o(n) is the Fokkcr-Planck operator in absence of flow (in equilibrium). In eq (|38l) . we have assumed that 
the shear flow direction and the shear gradient direction are x and y, respectively, and 7 is the shear strain which is 
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assumed to be sufficiently small. The distribution functions can be expanded into the power series of 7 and the Oil 2 } 
terms can be dropped when we calculate the shear relaxation modulus (because it is a linear response function[35|). 
From eqs (fT5|) and (|38|) . the distribution functions at t > can be calculated as follows. 



$(r,t) = $ cq (r) - e t£ °(°)^- [ 7 r y $c q (r)] 



+ rV £ ° (0) 

Kb I 



dv(r) 



dr x 

The dynamic equation for P(0,t) (eq (|14p ) can be rewritten as follows. 



dP(0,t) 
di 



i[i_P( , *)]_!. 

TO T 



P(0,<) 



(39) 



(40) 



The effect of the step shear on P(0, i) appears in the second term in the right hand side of eq (|40| only through 
<f>(r,t). However, from the symmetry, the integral over $(r,f) becomes 



/2B* 



,tC o {0) 



J dre- 3r2 / 2B ?»S>(r,t) 

= J dre- 3r2 ^ cq (r) + ^ J dre~^ 

= J dre- 3ra / 2 ^$eq(r) 
This means that a small deformation does not affect P(0,t). 

P(0,i) = P eq (0) 

From eqs (|13p. (j3"9")l and (|42p . we have the following expression for P(l,r,t). 



dv(r) 

dr x 



(41) 



P(l,r,t) = Peq(l,r) 



T_ e t[Co(l)-l/ro] 



R o 



• y 



^ P cq(0) [ fa' e (t-t')[C (l)-l/T ] 



k B Tr Q 



dr x 

e -3r 2 /2Rl e t'£a(0) 



dv(r) 
dr x 



^cq(r) 



(42) 



(43) 



The shear relaxation modulus is calculated from the average stress tensor as G(f) = a xy (t)/^. Substituting eq (|43p 
into eq (|35|) . we finally have the expression of the shear relaxation modulus G{t). 



G{t) = 



nopo 



R 



a t[£ (l)-l/r ] 



3k B Tr x | fo(r) 
^0 



dr.. 



r y P eq (l,r) 



£ p cq(°) [ dt > e (t-t')ic (i)-i/r ] 



TO 



-3r 2 /2Rl t'£ (0) 



dv(r) 
dr x 



r y $eq(r) 



(44) 



Intuitively, this result can be interpreted as follows. The first term represents the relaxation by the thermal motion 
and reconstruction of the network (bridges) . The second term represents the coupling between the spatial correlation 
of micellar cores and the bridge construction. That is, if the structure is deformed by the applied shear deformation, 
newly constructed bridge chains are also deformed. Such a coupling gives non-trivial relaxation modes. The second 
term in the right hand side of eq (|44p is missing in most of previous models. Eq f|44[) reduces to a single Maxwellian 
form in the limit of £ — > and £0 —> 0. This limit is, however, physically unreasonable because the network is no 
longer dense: The functionality becomes very small for £ — > and the network becomes sparse. Also, the correlation 
between micellar cores is not small and thus Co is not negligible. 

For a sparse network, the condition considered above turn to be reasonable. Generally the effective activity (fugacity) 
£ is small, and the Fokker-Planck operator in absence of flow is negligible compared with other contributions (A)(1) — ► 
0). The characteristic network size becomes larger than the characteristic interaction range between micellar cores, 
and thus the effects of spatial correlation between micelles become small. Therefore, here we consider the shear 
relaxation modulus of the sparse network model for the cases where ( < 1, A)(l) ~ 0) an d V ( r )/ K BT <C 1. The 
effective potential u(r) and the Fokker-Planck operator simply become 



u(r) 



3k B T 
2mRf 
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C(l,t)P(r)tt- — [ irx 5(t)P(r)] 



(45) 
(46) 
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and the shear relaxation modulus and characteristic modulus have simple forms 

G(t) « G e~ mt/T0 

dv(r 



dr 



3r x 



dr x 



^ cq (l,r) 



(47) 
(48) 



This result seems to be in harmony with the experimentally observed single Maxwellian type relaxation behavior 
of the telechelic polymer solutions with intermediate concentrations, in which sparse networks of superbridges are 
expected. 



2. Shear Viscosity and First Normal Stress Coefficient 



Under fast flows, telechelic polymer solutions show some nonlinear rheological properties. Here we consider the 
shear viscosity and first normal stress coefficient under shear flow. For the dense network model, they arc defined as 



r)(j) = ~ J dra xy (l)P ss (l,r) 



(49) 
(50) 



where 7 is the shear rate, and the shear flow and shear gradient directions are set to the x and y directions, as before. 
The quantities with the subscript "ss" represent the steady state quantities under shear. Experimentally, both 77(7) 
and ^1(7) are known to exhibit nonlinear behavior. (The nonlinear behavior of 77(7) and ^1(7) depend on various 
factors such as the polymer concentration}^ .) 

The steady state miccllar core distribution function is given as the solution of the following equation. 



£ ss (0)$ ss (r) = 



(51) 



We cannot obtain the explicit form of $ ss (r) because we have not specified the Fokker-Planck operator. (Even if the 
explicit form of £ S s(0) is given, generally it is still difficult to obtain the explicit form of $ ss (r).) However, some 
qualitative properties of $ S s(f*) can be argued without such detailed information. An important and general property 
of < & ss (t') is that it is anisotropic. This is because the operator £ S s(0) contains the advection term, which is intrinsically 
anisotropic. If the mobility tensor becomes anisotropic, the anisotropy will be enhanced. 

The steady state probability distributions P ss (0) and P ss (l,r) can be formally expressed in terms of $ ss (r). 



P ss (0) 



1 + e/ dre- 3r2 / 2R2 °$ BS (r) 



(52) 



P ss (l,r) = P ss (0)^T fe £ s fc s (l)[e^ 2 / 2ii o$ ss (r)] (53) 

fc=0 

Thus P ss (l,r) is also anisotropic under shear. Moreover, the operator £ S s(l) modulates the anisotropy in a rather 
complicated way. Therefore, the effect of the shear flow on P ss (l, r) is not trivial even in this simple case. 
From eqs (|4"9")h (|5P|) and (|53p . the shear viscosity and the first normal stress coefficient become 

riin) = ^f^Pss(O) / dr^±r^i(l)[e-^^ (54) 

Mi) - ^^PssiO) I *^-^f;^4 s (l)[ e -^/ 2 H^ ss(r)] (55) 

In general, both $ ss (t*) and £ S s(l) are anisotropic, and thus the dependence of 77(7) and ^1(7) on 7 is not simple. 

For the sparse network model, we can approximate the Fokker-Planck operator as C(l,t) « —jr y (d/dr x ). Then 
we have the following expressions for the shear viscosity and first normal stress coefficient. 

(56) 



-Pas(O) 



1+iJ dre- 3r2 / 2mR2 °$ l ™ ] (r) 
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, . s n cff p £k B T = /" 
7?(7) ~ ; ^ ss (0) / 



tq 7 



e -3r 2 /2mfl^[H (T , ) 



(57) 



#l(<V) 



n c sp £,k B T 



r 



+ 2 



P ss (0) / dr 



To 7 



tuRq 
( Z°l 



(58) 



Eqs (|57p and ([55)) may seem to be somehow similar to eqs ([Mf and ([53]) , and one may consider that the nonlinearities 
in the sparse network model are rather weak. Indeed, the origin of the nonlinearities itself is qualitatively the same 
as in the case of the dense network model. (The nonlinearities mainly come from the nonlinear and anisotropic 

dependence of & ss (r) on shear flow.) However, we should recall that 3>l!T'(r) is expressed in terms of convolutions of 
$ ss (r) (eq (|28p). and thus its nonlinearity will be enhanced. Therefore, the shear rate dependence deduced from eqs 
([571) and (|58p can be stronger than that in the dense network model. 



Here, we examine the shear thickening and thinning behavior. We expand the function 3?ss (r) into a Hermite 
polynomial series, following the previous workfl8j. 



e- 3r2 / 2mR o<Z>t ] (r) 



3/2 



-3r 2 /2mR l 



1 E ^,fc(7)2- 



(i+j+k)/2 



x Ht 



/2mR 



H, 



i,j,k=0 

\Z2mRo 



'2mR n 



(59) 



Here H(x) is the i-th order Hermite polynomial [33], and {Bi,j,k(j)} are the expansion coefficients. Substituting eq 
([59| into eqs ([56 ]) - ([58[) . we find that the shear viscosity and first normal stress coefficient are expressed in terms of 
{B iJ:k (j)} as 



77(7) 



7[1 + £Bo,o,o(7)] 



D 



1,1,0 



(7) 



— ) [-80,0,0(7) 
to 



2B 



0,2,0 



(7)] 



(60) 



*i(7) 



2n cS p £k B T 
7 2 [1 + ^0,0,0(7)] 

#1,1,0(7) 



Tq7 
m 



£2,0,0(7) - £0,2,0(7) 



— ) [£0,0,0(7) 
m 



2£o,2,o(7)] 



(61) 



From eqs (|BU| and ([6"Tjl . it is clear that the shear thickening and/or thinning are determined by the shear rate depen- 
dence of four coefficients, -Bo, 0,0(7), £1,1,0(7), £2,0,0(7), and -80,2,0(7)- vii) an d ^1(7) depend on these coefficients 
differently. Also, the 7-dependence of these coefficients is generally not simple. As a result of the competition be- 
tween these coefficients, several different types of nonlinear behavior can be reproduced. Because the anisotropy is 
enhanced in the case of the sparse network model, we expect that the sparse network model can exhibit variety of 
shear thickening and/or thinning behavior. (We can further analyze nonlinear behavior by expanding {-6^-^(7)} into 
power series of 7(3, although we do not show details here.) 



IV. DISCUSSION 



A. Shear Relaxation Modulus 



The shear relaxation modulus of the dense network model (eq (011)) does not reduce to the single Maxwellian 
relaxation as shown in Sec. IIII CTl As we mentioned, the conditions for the single Maxwellian relaxation (g cq {r) ~ 1 
and ^< 1 (or ^ < 1)) are satisfied if the concentration of micellar cores is low and the network is sparse. Experimental 
data shown in this work have suggested that the crossover concentration of our HEUR solutions is about c c ~ 4wt% 
(cf. Figures [5] and [5]) . For c c < 4wt% we can employ the sparse network model and thus the single Maxwellian 
relaxation is reproduced. For c c > 4wt%, we should employ the dense network model. 
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The expression of the shear relaxation modulus in the dense network model is not so simple. To make the expression 
analytically tractable, we replace the Fokker-Planck operators in eq by constants, as £o(n) — > — 1/f (n). Such 
a replacement corresponds to an approximation for the operators by their characteristic eigenvalues. With this 
approximate replacement, eq (|44[) reduces to a simple form. 

, , e -t/f(0) _ g-i[l/f(l) + l/T ] 

G(t) « G e-*[^(D+ Vro] + Gl e e (62) 

1 + t /t(1) -r /r(0) 

Gi is a sort of characteristic modulus defined as follows. 

n o p o ^P cq (0) f 3r x r 2 y dv{r) 3 r 2 /2^, 



ft, ^i*-rv ^ (63) 

(Gi is expected to be small compared with Go-) The first term in the right hand side of eq (|62p is the single Maxwell 
type relaxation. The characteristic relaxation time is given as tqt(1) /[t® + r(l)]. (This relaxation time is longer 
than the intrinsic time To and changes with the polymer concentration.) The second term in the right hand side of 
eq (|62p represents to the deviation from a single Maxwellian form, and it originates from the coupling between the 
spatial distribution of miccllar cores and the bridge construction. The characteristic relaxation times f (0) and f (1) 
are expected to be larger than To- In the short and long time limit, the deviation term asymptotically becomes 

t/f (0) _ p -t[l/f (l)+l/ro] (t/T + O((t/T ) 2 ) (*<r ) 

e -*/ f (°) ^ (64) 



l + To /f(l)-ro/f(0) — — — - (t»r ) 

I 1 + t /t(1) - t /t(0) 

Intuitively, the deviation term slightly decelerates the relaxation in the short time scale and add a slow and weak 
relaxation mode in the long time scale. As a result, the relaxation spectrum becomes broader than the single 
Maxwellian spectrum. 

As an example, we show the storage and loss moduli, G'(cj) and G"(w) calculated from eq (|62)l for some parameter 
values in Figure The deviation from the single Maxwellian behavior is clearly observed. Similar trend can be 
observed in the experimental data, Figure We expect that this is because our dense network model captures 
the essential physics qualitatively. However, we should notice that the estimate shown here is rough and cannot be 
quantitatively compared with experimental data. In reality, the Fokker-Planck operators give multiple relaxation 
times and the shear relaxation modulus (or the storage and loss moduli) will be much broader. Here it may be 
worth mentioning that the deviation from the single Maxwellian is already predicted in the transient network model 
by Tanaka and Koga[38j|. However, in their model, the storage and loss moduli data become sharper (in the mode 
distribution), whereas our model and the experimental data show the broadening. This difference could result from 
the lack of the correlation between micellar cores in the Tanaka-Koga model. 



B. Concentration Dependence of Characteristic Time and Modulus 



Experimental data show clear concentration dependence of the characteristic relaxation time r and modulus Go- 
In the single chain transient network models without spatial correlation effect and the network functionality, the 
experimentally observed concentration dependence cannot be reproduced. 

We consider the dependence of the characteristic relaxation time r in the conventional single chain models. The 
characteristic relaxation time r is just the same as the average bridge destruction time to, t = To. to is considered 
to be determined by the characteristic relaxation time of a polymer chain (such as the Rouse and Zimm times) and 
the activation energy (energy barrier). Naively, both of them are almost independent of the polymer concentration 
c, and we can roughly estimate as t cx c . This estimate for the concentration dependence of r is not consistent with 
experimental data. (Although it may be possible to remedy the model by introducing some concentration dependence 
to parameters such as To, such an ad-hoc approach is not always physically reasonable.) 

In our models, the concentration dependence of t is not that simple. For the dense network model, due to the cou- 
pling between the association-dissociation dynamics and the dynamics of miccllar cores, the characteristic relaxation 
time has been estimated to be t w tot(1)/[to + t(1)] (see cq (jMI)- This relaxation time can depend on c through 
f(l). For the sparse network model, from eq (|47p the characteristic relaxation time is given as t = To/m, which can 
also depend on c because m is a function of c. Besides, the c-dependence of r of the dense and sparse network models 
is expected to be different (because of the difference between the c-dcpcndence of f (l) and to). This is qualitatively 
consistent with the experimental results (Figure [S^a)). 
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C. Fast Mode 

The storage and loss moduli data of the HEUR solutions exhibit the fast mode at the high frequency region. Our 
models do not take into account the local relaxation dynamics of polymer chains and thus they cannot reproduce the 
fast mode. Nonetheless, we can estimate how the fast mode behaves for the dense and spare network cases. 

In micellar network structures, chain ends of associative polymers are connected each other and thus the relaxation 
mode distributions are given as the eigenmode distributions of networks. This situation is similar to the Rouse model, 
where the linear viscoelasticity is described by eigenmodes of linearly connected springs. For the sparse network, we 
expect that the fast mode mainly reflects the relaxation of a superbridge structure. In analogy to the Rouse model, the 
relaxation time depends on the number of chains per superbridge m as Tf as t oc to 2 tpeo- Here tpeo is the relaxation 
time of a free PEO chain. Thus the relaxation time of the fast mode would be longer than the relaxation time of 
the PEO solution. Moreover, the viscosity depends on m as f?f as t oc v c ^m? oc m, with v Q s oc mT x being the effective 
number density of supcrbridges. Thus G"(u>) of the HEUR solution is expected to be larger than the PEO solution. 
Although the estimate above is very rough, the trend is consistent with experimental data (Figure [T]) . 

For the dense network, the situation is qualitatively different. In the dense network, polymer chains are connected 
in a rather complicated way. The eigenmode distributions will be much broader than the case of the sparse network, 
and the relaxation time will be much longer. As a result, the fast mode of the dense network will be broader than 
the sparse network. Also, the relaxation time of the fast mode will be longer than one of the sparse network. If we 
assume that networks are fractal structures, the storage and loss moduli are described by the power law (just as in 
cases of gels[3^|). The power law exponent is expected to decrease as the concentration increases and the network 
becomes denser (well-percolated) . This is consistent with experimental data shown in Figure 2] Therefore we can 
again conclude that, for both the spare and dense networks, the fast mode reflects the relaxation behavior of main 
chains in the connected (associated) network structures. 

D. Nonlinear Elasticity and Stretch Dependent Reconstruction Rate 

Here we discuss the effects of the nonlinear elasticity and the stretch dependent reconstruction rates in detail. In 
particular, wc focus on mechanisms for the shear thickening and thinning, for which these factors have been assumed 
to be important. 

In some of conventionally examined transient network models, the shear thickening and thinning are explained as a 
result of the competition between two factors; the nonlinear elasticity and the stretch dependent dissociation ratefliT- 
lvi\ . Under fast flow, the nonlinear elasticity increases the shear stress whereas the stretch dependent dissociation 
decreases the number of highly stretched chains and thus decreases the shear stress. As we have shown, however, even 
without the nonlinear elasticity and the stretch dependent reconstruction rates, some nonlinear rhcological behavior 
can be naturally reproduced. In our models, the shear stress is strongly affected by the steady state distribution 
function 3> ss (f*) given by eq (|51[) . Thus, the structural anisotropy can be the main mechanism giving rise to the shear 
thickening and/or thinning in actual telechelic polymer solutions. 

This implies that the nonlinear elasticity and the stretch dependent dissociation rate are not always essential 
ingredients of transient network type models. There are many different possible dynamic models which can reproduce 
the same shear thickening and/or thinning behavior. Even if a specific model reproduces the shear thickening and/or 
thinning behavior of a specific sample perfectly (with some parameter fittings), it does not guarantee that the model 
is physically correct. Actually, a recent experimental work showed that the effect of nonlinear elasticity is negligibly 
small for typical HEUR solutions, based on a simple energy balance argument [l8j]. Careful and systematic comparison 
of theoretical predictions and experimental data are absolutely necessary for validating a theoretical model. 

Some conventional models predict that the average bridge chain is smaller, in spatial size, compared with the 
average loop chain size and/or dangling chain size in equilibrium. (The fact that the stretch dependent dissociation 
rate affects the equilibrium chain statistics has already pointed by Tanaka and Edwards]?}-) However, judging from 
molecular simulation data[22[ or scattering data 0, [lfj, the average bridge chain size is larger than the loop/dangling 
chain size. Moreover, if the average bridge chain size is much smaller than the loop or dangling chain size, uniformly 
spanned networks cannot be formed. Under such conditions, some simple stretch dependent dissociation models will 
become physically odd. This consideration suggests that the spatial correlation (spatial structure) is an important 
factor determining the rheological behavior of telechelic polymer solutions. Of course, the stretch dependence may 
be important for some cases. The reassociation process in telechelic polymer solutions needs to be modeled with very 
careful considerations. 

It is fair to mention that our models still lack various features which may be important for some cases. For 
example, the direct bridge-to-bridge transition is not included in dynamic equations (|13[) or (|3ip . The bridge-to- 
bridgc transition can contribute to several rheological properties, especially when a network is dense and the number 
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density of micellar cores is high. The higher order spatial correlation functions, being important for such dense 
systems, are not explicitly considered in our models. If we rigorously construct the dynamic equation for the two 
body correlation in dense systems, generally we have the three body correlation function. (This is qualitatively the 
same as so-called the BBGKY hierarchy in the liquid state theory 3l|.) Also, the interactions between bridge chains 
or micellar cores are considered in our models only via the mean field type approximation, which may not be always 
acceptable. The polydispersity of molecular weights or aggregation number are also ignored in our models. These 
features are beyond the scope of this work and are not further discussed here. Elaborated theoretical investigations as 
well as detailed molecular simulations will be required to understand the molecular mechanisms in telechelic polymer 
solutions. 



V. CONCLUSIONS 



We conducted linear viscoelastic measurements for HEUR solutions with different polymer concentrations. We have 
shown that the HEUR solutions exhibit fast and slow relaxations, and both of them depend on the concentration. 
Our HEUR solutions exhibited different rheological behavior at low and high concentrations. These results suggest 
existence of two different concentration regimes for telechelic polymer solutions. 

To explain these experimental results, we have considered single chain type transient network models for dense and 
sparse networks, which take into account the information of spatial correlations between micellar cores and the network 
functionality. The spatial correlation largely affects the bridge construction dynamics in our models. We considered 
the linear and nonlinear rheological properties of the models, in absence of the conventionally considered mechanisms, 
such as nonlinear elasticity and the stretch dependent bridge destruction rate. The sparse network model reproduces 
the well-known single Maxwcllian type shear relaxation modulus for solutions of low polymer concentrations. On the 
other hand, the dense network model gives non-single Maxwcllian type shear relaxation modulus, which is consistent 
with experimental results for high concentration HEUR solutions. 

For the steady shear viscosity and first normal stress coefficients, both the dense and sparse network models 
exhibit nonlinear behavior such as the shear thickening and thinning. The nonlinear rheological behavior originates 
mainly from by the steady state distribution function of micellar cores, $ ss (r). This is consistent with a recently 
proposed anisotropic bridge formation model [lU- Our models suggest that the nonlinear rheological properties are 
also strongly affected by the polymer concentration. The spatial correlation is one possible origin of the nonlinear 
rheological behavior. The measurements and analyses of nonlinear rheological properties for the HEUR solutions are 
now in progress, and will be reported elsewhere [36j. 
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Appendix A: Anisotropic Bridge Formation Model 



In this appendix, we briefly show the anisotropic bridge formation model proposed in Rcf (18j |. We also show that 
the sparse network model formulated in this work reduces to the anisotropic bridge formation model under some 
conditions. 

First we briefly show the anisotropic bridge formation modelfl8|. In this model, the number fraction of bridges 
(superbridges) is assumed to be constant even under shear. We introduce the normalized probability distribution for 
an end-to-end vector, ip(r,t). The dynamic equation under a steady shear flow is described by a birth-death type 
master equation 
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^ = -ir y d -^l + I*(r,t) - U(r,t) (Al) 
at dr x r r 

where we have set the shear flow and shear gradient directions to the x and y directions, respectively. 7 is the shear 
rate, fo is the characteristic bridge reconstruction time, and (f>(r, t) is the distribution function for a newly constructed 
bridge (superbridge) . <f>(r,t) can be interpreted as a source function. In general, <j)(r,t) depends on the shear rate 
7, and the nonlinear rheological properties are mainly determined by <fi(r,t). The detailed dynamics of <fi(r,t) does 
not need to be fully specified. (In Ref [l8[, only some expansion coefficients at the steady state are utilized for the 
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analysis.) If the clastic potential of the end-to-end vector is given as a harmonic form, the stress tensor can be 
expressed simply as 



<j = vksT I dr 



3rr 



ijj(r) 



(A2) 



where v is the number density of elastically active bridges and ro is their average size. From eqs (|A1[) and (|A2[) . one 
can analyze some rheological properties. 

Now, we show that our sparse network model can reduce to the anisotropic bridge formation model. Because the 
polymer concentration of the HEUR solution used in Ref. fhsj is lwt%, the experimental system is in the sparse 
network regime as judged from the crossover concentration in Figures [3] and [5] (c c sa 4wt%). We define the normalized 
bond vector distribution function as 



if>(r,t) 



P(l,r,t) 
1 - P cq (0) 



and the normalized source function as 

6(r,t) = - ) , ' 
^ ' ' W(0\l,r) 



](r,t) 



1 - Pcq(O) 



1 - Peq(O) 



(A3) 



(A4) 



Eq (|A4[) means that the source function is determined from the distribution function of miccllar cores, $(r,£). (At 
the steady state under flow, the source function (pssif) is determined by § ss (r) in eq (|5ip. Namely, the steady state 
Fokkcr-Planck operator C ss (0) determines 4> ss (r).) Substituting eqs (| A3|) and (|A4[) into eq (f3~Tj) . we have the dynamic 
equation for ip(r,t). 



^M=£(l l t)^(r,t) + l^(r,t) 
at r 



TQ 



-ip(r,t) 



(A5) 



where we have set To = tq/iti. Further, if we neglect the Brownian motion and the interaction between micelles and 
set C(l, t)ip(r) = —jr y d^(r)/dr x , eq (|A5[) reduces to cq (|A1[) . Also, the stress tensor is given as 



tr = n e ffPo[l - -Peq(O)] / dr 



du(r) 
dr 



r - k B Tl 



ijj[r) 



(A6) 



With the harmonic form of the effective potential (cq (|45|) with ro = s/mRo), cq (|A6p becomes essentially the same as 
eq (|A2[) . Thus the sparse network model formulated in this work reduces to the anisotropic bridge formation model. 

However, it should be noticed that the dynamics of the source function is not explicitly specified in the anisotropic 
bridge formation model [l^. This means that we can use eq (|A1|) even when the source function is not given by eq 
141). There are many other possible molecular models which reduce to the anisotropic bridge formation model. 
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Figure Captions 

FigureQ] Master curves of storage and loss moduli for the lwt% HEUR aqueous solution. The reference temperature 
is T r — 25° C and the superposition is performed for the slow (low frequency) mode. 

Figure [2] Master curves of storage and loss moduli for the lwt% HEUR aqueous solution. Data are the same as in 
Figured] but shifted with determined for the lwt% PEO solution. The reference temperature is T r = 25°C. For 
comparison, the loss moduli for the lwt% PEO solution at T = 25°C is also plotted. 

Figure [3] Horizontal shift factors of the slow mode of the HEUR solution and the fast mode of the PEO solution, 
ar,o and ax,i- The reference temperature is T r = 25°C. 

FigureUJ Concentration dependence of storage and loss moduli for the HEUR aqueous solutions at T = 25°C. The 
concentrations are c = 1, 2, 5, and 10wt%. 

Figure Storage and loss moduli for the HEUR aqueous solutions with various concentrations at T = 25° C. 
Angular frequency and moduli are rescaled by the inverse characteristic time 1/r and the characteristic modulus Go 
for the slow mode, respectively. 

Figure [6j Concentration dependence of (a) the characteristic time r and (b) the characteristic modulus Go at 
T = 25°C obtained for the HEUR aqueous solutions. Broken and solid lines represent the fitting results for low and 
high concentrations, respectively. 

Figure [7] A schematic image of the mean field single chain model for a dense network (the dense network model), 
(a) a dense network formed by elastically active chains (solid black curves) and micellar cores (circles), (b) a local 
structure in a dense network. Solid and dotted curves represent elastically active and inactive polymer chains, 
respectively. Circles represent micellar cores. The dense network model is formulated for a single tagged chain in the 
system. There is the effective interaction between micelles, which is expressed by the potential v(r). A bridge chain 
feels the elastic potential u(r). 

Figure |5J A schematic image of the mean field single chain model for a sparse network (the sparse network model), 
(a) a sparse network formed by elastically active chains (solid black curves) and micellar cores (circles), (b) a local 
structure in a sparse network. Solid and dotted curves represent elastically active and inactive polymer chains, 
respectively. Circles represent micellar cores. The sparse network model is formulated for a single tagged supcrbridgc 
in the system. v(r) and u(r) are the effective interaction potential between micelles and the elastic potential of a 
supcrbridge. 

Figure M Storage and loss moduli (G'(w) and G"{uj)) deduced from a dense network model for f = f(0) = f(l) 
calculated under some approximations. For comparison, G'(cj) and G"(cj) of a single Maxwellian model (Gi/Go = 0) 
are also shown. 
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